The spectral edge of some random band 

matrices 

Sasha Sodin* 
June 25, 2010 

Abstract 

We study the asymptotic distribution of the eigenvalues of random 
Hermitian periodic band matrices, focusing on the spectral edges. The 
eigenvalues close to the edges converge in distribution to the Airy point 
process if (and only if) the band is sufficiently wide {W ^ N^^^). 
Otherwise, a different limiting distribution appears. 

1 Introduction 

In this paper, we study the edge of the spectrum of random Hermitian pe- 
riodic band matrices. The N x N Hermitian random matrix H = Hn with 
rows and columns labelled by elements of li/NX has independent entries 
above the main diagonal, and 

Huv = if \u — v\j\f = m.m{\u — v\, N —\u — v\) > Wf^ or u = v. (1.1) 

To simplify the exposition, we assume that either 

P{i7„, = 1} = P{i/,, = -1} = 1/2 , {]<\u-v\m<Wn (1.2) 
( "random signs" ) , or 

Huv = exp{iUuv) , Uuv ~ U{0, 27r) , < \u - v\n < Wn (1.3) 



*Tel Aviv University, sodinalc@post.tau.ac.il. Supported in part by the Adams Fellow- 
ship Program of the Israel Academy of Sciences and Humanities and by the ISF. 



1 



("random phases"), and defer the discussion of possible generahsations to 
the last section. For the same reason, we assume that W = Wn — )■ oo as 
N ^ oo. 



The matrix is closely related to the graph S = {'E/N'E, £), where 

{u, v) 0<\u-v\<Wn , (1.4) 

and can be viewed as a Hamiltonian of quantum evolution in a disordered 
environment on 9 (f ll.Sp corresponds to broken time- reversal symmetry.) We 
refer the reader to the work of Fyodorov and Mirlin [TT] for a thorough 
discussion of physical motivation. 



Random matrices similar to Hn have been studied in both mathematical and 
physical literature. We survey some of the results that pertain to the current 
work, and refer to [H] and to the recent review of Spencer pB] for detailed 
bibliography. 

(A) Bogachev, Molchanov, and Pastur [1] have proved that the empirical 
spectral measure (or integrated density of states) 

Ativ(a) = # jeigenvalues of ^^y- "]} ^^-^"^ 

converges (weakly, in distribution, as — )■ oo and Wj\f — >■ oo), to the (deter- 
ministic) Wigner measure 

/a 2 
-(1 - aj)]^'^dai . (1.6) 
-oo 

That is, the global behaviour of the spectrum of if at is similar to that of 
Wigner matrices (which correspond to the special case Wn = N/2.) 

(B) One of the interesting questions concerning local eigenvalue statistics is 
the crossover between the Random Matrix regime and the Poisson regime. 
The Thouless criterion [2S], applied to random band matrices by Fyodorov 
and Mirlin [11], predicts the following. 

If the mixing exponent pmix of the (classical) random walk on S is much 
larger than the eigenvalue spacing near ao: 

mean spacing at ao 

Pmix ^ ~, TT 7 5 

mean density at 
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then the eigenvalues of near obey Random Matrix statistics, and the 
corresponding eigenvectors are extended. If 

mean spacing at oq 

Pmix ^ TT 7 ) 

mean density at ao 

the eigenvalues near ao obey Poisson statistics, and the corresponding eigen- 
vectors are localised. 

For our graph S, Pmix is of order W^/N'^, and in the bulk of the spectrum, 
— 1 < «o < 1; (11. 6p suggests that 

mean spacing at ao ^ ^^j^ 
mean density at ao 

Thus the crossover should occur at x \/N. 

On the physical level of rigour, these predictions have been justified by Fyo- 
dorov and Mirlin [TT], who have derived a detailed description of both asymp- 
totic regimes (and also of the crossover.) So far, these results resist math- 
ematical justification (see however the work of Khorunzhiy and Kirsch [16] 
for some relatively recent developments.) 

We remark that Spencer and Wang [331 Ch. Ill] have managed to make 
one direction of (a slightly different form of) the Thouless criterion rigor- 
ous under certain assumptions, in a fairly general setting. However, for now 
these assumptions have not been verified for the problem under considera- 
tion. We refer the reader to the review of Spencer |28j for a discussion of the 
mathematical approach to the Thouless criterion and its application to band 
matrices. 

(C) We focus on the edge of the spectrum: ao = ±1. There, 

mean spacing at ao ^ N^'^^'^ N~^^^ 
mean density at ao A^~^/^ ' 

therefore the crossover should occur at x N^^^. On the physical level of 
rigour, this has been confirmed by Silvestrov [23]. In fact, one of the aims of 
the current paper is to put some of the methods and results of [23] on firm 
mathematical ground. 

Now we state the main results. 
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© For Wn > A^s/e (i.e. Wn/N^'"^ — > +00), set 

= # jeigenvalues of ^^^= in (l - +00] | , 

a.(A) = # {eigenvalues of in ( - 00, -1 + ^] } . 

Theorem 1.1. IfW^/N^^^ 00 as N ^ 00, the measures a r{X) andaL{\) 
converge in distributioi^ to 2ti/3(— A), where 211^ is the (random) distribution 
function corresponding to the Airy point process, /3 = 1 for and (3 = 2 

for m. 

Hence, according the Tracy-Widom theorem IS^ lSl ^ . the scaled extreme 
eigenvalues 2N'^l^{a^^-l) and -2N'^/^{a^^,, + l) of Hn/ {2y/2W^) converge 
in distribution to the Tracy-Widom law TW^. 

2) For K < N^/^, set 
o'r(X) = — # < eigenvalues of — 7^= in f 1 -77, +00I 



= —77— # I eigenvalues of — in ( — cx), — 1 H 

AT ^1 ^ 2V2W^ ^ 2Wf^ 

Theorem 1.2. If Wn — ?• 00 and W^/N^^^ 0, the measures and 
converge in distribution to the (same) deterministic measure ajs. We have: 

o-^(A) = |-A3/2 + 0(A) , A^+00 , 
ovr 

cTp{X) < C exp{-C\Xf^) , A ^-00. 

In some sense, the eigenvalues close to the edges behave as independent ran- 
dom samples from o"^. The method of the current paper yields the following 
manifestation of this belief: 

P {2Wf{l - > a} - exp i -^a^(A) I = o(l) 



-'^on test functions ./ G C(Il) such that supp/ n R+ is compact 
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uniformly in A G R. Unfortunately, our description of the left tail of ap 
is not sufficiently precise to deduce convergence to a max-stable law (cf. 
Gnedenko \T2\ for the description of domains of attraction of max-stable 
laws.) 

Khorunzhiy [15] has proved (for a slightly different class of band matrices) 



as — )■ oo (actually, he has established a stronger form of convergence.) 
He has conjectured (private communication) that the same conclusion holds 
under the weaker assumption Wn ^ logA^. We confirm this conjecture: 

Theorem 1.3. IfWN/logN +00, then \\Hn / i2^/2mv)\\ 1. 

As one can see from the argument of Bogachev, Molchanov, and Pastur 
this result is sharp, meaning that the conclusion fails if Wat/ log — )■ 0. 

The proofs of the three results are based on a modification of the moment 
method. The moment method has been applied to random matrices since 
the work of Wigner pTj. Bogachev, Molchanov, and Pastur |1] have applied 
it to study the spectrum of band matrices (see above.) 

The moment method appears particularly useful to study the eigenvalue 
statistics at the spectral edge. In [26], Soshnikov has applied it to derive 
the limiting distribution of the extreme eigenvalues of Wigner-type random 
matrices. Extensions of his approach have allowed to solve a number of 
related problems, e.g. [IHIEIIISZ] • 

In this paper, we apply a modification of the moment method, which 
goes back at least to the work of Bai and Yin [2] (see [25] for more detailed 
references.) In [8], it has been used, in particular, to give another proof of 
Soshnikov' s result. The combinatorial technique of [8] is the main ingredient 
of the current work. In Section [2l we review this technique, and formulate 
the combinatorial statements needed to prove the main results. 

Another ingredient of the proof is an asymptotic description of the (classical) 
random walk on 9- We prove the necessary facts in Section [31 In Section HJ we 
apply these facts to count subgraphs of S of a certain form. In Section we 
specialise to the setting of Section |2l and prove the combinatorial statements 
formulated there. 



that, if Wn > log^/2 N, 




D 
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(1.7) 
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Section E] collects some facts related to Levitan's uniqueness theorem [T7] 
which we use in the sequel. 

In Section [TJ we conclude the proofs of the main results. Section |H] explains 
how to modify the proofs written for fll.2l) to deal with fll.3p . In Section [HI 
we discuss generalisations and related open problems. 

Notation: In this paper, C, c, • • • stand for positive constants, the value of 
which may change from line to line. If X, Y are quantities depending on some 
large parameter, X < F <^ F > X X = o(Y) X/Y 0. 



2 Preliminaries 

From this section, we omit the subscript N, and write H for and W for 
Wn- Also, we consider the matrices with entries (11. 2p : in Section |8] we shall 
explain the modifications needed for the case fll.3p . 

Let Un{cos6) = sin((n + 1)0)/ sm6 be the Chebyshev polynomials of the 
second kind; set f/_2 = f/-i = 0. Let 

^ ' \ \2^2W-l) 2W-1 \2^2W-l) 

The following lemma (see e.g. [25] for a more general version) is at the basis 
of our considerations. 

Lemma 2.1. For any Hermitian N x N matrix H satisfying 



±1, Q<\u-v\n<W 
0, otherwise 



and any uo,Un € 'Z/N'Z, 

where the sum is over all {n — l)-tuples {ui,U2, ■ ■ ■ ,Mn-i) (which we regard 
as paths Pn = UqUi ■ ■ ■ Un), such that 

(a) < \uj — Mj+iIat < W , < j < n — 1 (pn is a path on S ), 

(b) Uj, < j < n — 2 (pn is non-backtracking.j 
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The following corollary is immediate from the lemma and fll.2p . 

Corollary 2.2. For the random matrix H as defined in the introduction, 

Etri7W^)Hri7W2))...^^^(n(fc)) 

is equal to the number of k-tuples of paths (shortly: k-paths) 

Pn{l),- ,n{k) = u\u\ ■ ■ ■ 'U^(^) ulul ■ ■ ■ ^^(2) ' ' ' ""o""! ' ' ' ^n{k) 

that satisfy (a),(h), and also 

(c) = Uq, ^ ^ j ^ k (the paths are closed); 

(d) the number 

# {{hj) Wi=u, = t;} - # {{i,2)\u\ = V, M^+i = u} 
is even, for any u,v & li/NX. 

As in [8j, we group the /c-paths satisfying (a)-(d) into topological equivalence 
classes, which are in one-to-one correspondence with k-diagrams: 

Definition 2.3. Let /3 G {1,2}. A k-diagram is an (undirected) multigraph 
G = (y, E), together with a fc-tuple of circuits 

p = ulu\---ul, ulul---ul, Uqu'1---Uq (2.1) 

on G, such that 

• p is non-backtracking (meaning that in every circuit no edge is followed 
by its reverse, unless the edge is a loop uu); 

• For every (u, v) G E, 

# {{hj) I u] = M, =v} + u] = V, = u} =2; 

• the degree of Uq in (5 is 1; the degrees of all the other vertices are equal 
to 3. 



7 



Remark 2.4. We emphasise that G is a multigraph in which the coinciding 
edges are distinguished. Thus, strictly speaking, a circuit is not uniquely 
determined by the vertices it passes. Still, we find it convenient to use the 
notation as in Ii2.1\) . Next, we do not distinguish two diagrams which are 
isomorphic in the natural sense. Thus, by a diagram we actually mean an 
equivalence class (e.g. in the second part of the following lemma.) 

The following lemma summarises some properties of /c-diagrams from |8l 
Part II]. 

Lemma 2.5. 

1. For every k-diagram, there exists an integer s > k ( "non-orientable 
genus", cf. Figures 1,2), such that the diagram has 2s vertices and 
3s — k edges. 

2. The number Dk{s) of k- diagrams corresponding to a given s satisfies 




Figure 1: Some 1-diagrams: s = 1 (left), s = 2 (center, right) 
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Thus our goal is to compute the number of /c-paths corresponding to a given 
diagram. It will be convenient to consider connected diagrams only. This 
can be done as follows: set 

T(n(l)) = Etr/J("(i)) , 

T{n{l),n{2)) ^ EtT H^^^'^hr H^""^^^^ - T{n{l))T{n{2)) , 



k 

T(n(l), ■ ■ ■ , n{k)) = EHi/W^)) -J^U ^^({^O') W) , 

j=i n Pen 

where the sum is over non-trivial partitions of {1, 2, ■ ■ ■ , k} into disjoint sets. 
It is not hard to see that indeed T(n(l), • • • , n{k)) counts the number of k- 
paths Pn(i), -- ,n{k) that satisfy (a)-(d) and correspond to connected diagrams. 

Now we formulate the two main technical statements of this paper. 
Proposition 2.6. Let W > N^/^, R>0. For any 

1 < n(l) < n(2) ■ ■ ■< n{k) < RN^^^ 
such that n{l) + ■ ■ ■ + n{k) — 2n is even, 

Tjnjl),--- ,n{k)) _ 
{2W - 1)" 

k 

4>u{n{l)/N''\ n{k)/N'/') + N'/h^inil), n(^)) , 



b=i 



where 



1. (pkizi, ■■■ ,Zk) = Yls>k9k,s{zi, ■■■ , Zk), gk,s being a continuous homoge- 
neous function of degree 2){s — k), 

gkA^ir-- ,zk) < ^ " 



2. eAr(n(l), ■ ■ • ,n{k)) = o{l), where the implicit constant depends only on 
R and on W/N^'^, and 

sr- £jv(n(l),--- ,n(/c)) 

^ n{l)---n{k) - 

l<n(l)<-<n(fc)<iVi/3 ^ ^ ^ ^ 



9 



If n{l) + ■ ■ ■ + n{k) = 2n + l is odd, T(n(l),--- ,n{k)) = 0. 
Proposition 2.7. Let W < N^^^ , R>0. For any 

1 <n(l) < ■■■ <n{k) <RW^/^ 
such that n{l) + ■ ■ ■ + n{k) = 2n is even, 
r(n(l),---,n(A;)) 



1^6/5 

where 



{2W - 1)" 

k 

b=i 



1. ipk{zi, ■■■ ,Zk) = Y.s>k hk,s{zi, ■■■ ,Zk), hk,s being a homogeneous func- 
tion of degree (5s — 5k — l)/2, continuous outside the origin. 



hk,s{zii 



{c\\z\\r 



, , 3s — 5fc-f-l 
[CS] 2 



2. ENijiil), ■ ■ ■ ,n{k)) = o(l), where the implicit constant depends only on 
R and on N^/^/W, and 



£Af(n(l), ■ ■ ■ ,n{k)) 
n{l) ■ ■ ■ n{k) 



l<n(l)<---<n(fc)<VK2/5 

If n{l) + ■ ■ ■ + n{k) = 2n + lis odd, T(n(l),--- ,n{k)) = 0. 



3 Random walk on the circle 



Let S = (^jNlj, £) be defined by (11.41) (in this section we do not assume that 
^ 1.) Denote by W„(i?) the number of paths of length n in S between 
two (fixed) vertices u,v such that \u — v\n = R. We prove the following 
statements. 



Proposition 3.1. // K n < N^/W^, R < n^/^W , and R < 0A9N, th 



en 



{if- (1+0(1)) [^(«' + l)(2«'+l) 



-1/2 



exp 



n{W + l){2W + l) 
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Proposition 3.2. IfN'^/W'^ < n, 

Proposition 3.3. Without any assumptions on n, N, R, W , 
W„(i?) 



{2wy 



< c 



For fixed W , Proposition 13.11 follows from Richter's local limit theorem for 
lattice variables (see Ibragimov and Linnik, [131 Ch. VII].) However, we 
need the asymptotics to be uniform in W , which we have not found in the 
literature. Proposition 13.21 (perhaps with an extra logarithmic factor in the 
assumptions) follows easily from the spectral estimates on the mixing time. 

Let A be the adjacency matrix of S- As S is invariant under cyclic shifts, 
the discrete Fourier transform diagonalises A. We state this as a lemma: 

Lemma 3.4. {2W)^^A = X]^o^ ^kfk ® fk, where 

sin(l^^) / ^nk 
ak= ffc cos iy + 1 — 



an 



dUi)=exp{^-^}. 



Proof. The vectors are of unit length, hence we only need to check that 
Afk = 2Wakfk- Let w = exp ^. Then 

{Afk){m)= Yl ^'' = him) 

0<\e-m\N<W 0<\£\n<W 

Now, 

,l_^{2W+l)k 



1 -w'^ 

o<\e\N<w 



UJ 



{W+l/2)k _ ^k/2 ^-k/2 _ ^-{W+l/2)k 



UJ^I"^ — UJ^^/'^ 
ljWk/2 _ ^^-Wk/2 



^k/2 _ ^-k,2 yuj^"^^'^"' + c.-(^+^)^-/2) = 2Wa, . 

□ 
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Denote by Sq, - ■ ■ , ^at-i the standard basis in 



. Then 



Sj = N ' exp — fk , 



A;=0 



therefore 



^ = {{A/2Wr5R, 5o) = iV-i 5^<exp 



27iiRk 



Thus we have to estimate the above sum. Informally, the argument is as 
follows: max(|ai|, ■ ■ ■ , |a7v-i|) = 1 — Qi^N"^ /W^), whereas Qq = 1. If n 3> 
N^/W"^, the sum is dominated by the first addend, which is equal to 1. 
If n <^ N'^/W'^, the sum can be replaced with an integral, which is then 
evaluated using the saddle point method. And now to the formal proof. 

Denote 

sin [Wttz] 



cos [{W + l)nz] . 



W sin [ttz] 

The following lemma summarises some properties of /: 
Lemma 3.5. The function f is an entire function with period 1; 

< exp {-cW'x'} , -1/2 < X < 1/2 . 

\f{W)\ 

Proof of Proposition \3.1i Choose ^^^^^ <^ r] <^ 1, and consider two cases, 
(a): R < r|^/nW. Then, according to (13.11) and Lemma [3.51 



{2wy' 



N 



1 + 



-N/2<k<N/2 
-N/2<k<N/2 

The contribution of the second addend is negligible: 



r]\k\y/nW\ 



-N/2<k<N/2 



■r]\k\y/nW 



N 



< 



N 



-1 



-N/2<k<N/2 



cW'^nk'^ \ r]\k\y/nW 



iV2 



N 
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Then, 



+ E ([/(^/^r - {-^(W^ + l)(2W^ + 1)^}) 
The second sum is neghgible, whereas the first one is 

(l + o(l))iV-i E exp|-^(l^ + l)(21^ + l)— I 

k=—oo ^ 



nn 



;i + o(i)) _(w^ + i)(2iy + i) 



-1/2 



according to the Poisson summation formula. 

(b): -R > rjy/nW . According to (13. ip . Lemma [3.5[ and the residue theorem, 



{2wy 



-N/2<k<N/2 



N 



^•'^ 1- exp[2mNz] 



where the contour C encloses the zeros k/N, —N/2 < k < N/2, of the 
denominator, as in Figure [3l 



y+ 



Figure 3: The contour. 



The sum of integrals along the vertical parts of C vanishes, hence 

Jiy+-l/2 ' Jiy.-l/2 
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The first integral: 

= / exp0+(2;)- , 

iy+~i/2 Jiy+~i/2 1 " exp [2mNz\ 

where (f)+{z) = nlnf{z) + 2T^iRz. We choose ?/+ to make (f)'^ vanish at 
Zj^ = iy^. Expanding (p^ in Taylor series at 0, we obtain: 

{W + 1){2W + V)z^ (1 + 0{z^W'^)) + 27rii?z . 



3 

This expansion is of course different iable, so 

(\^^{z) = —{W + 1){2W + l)z (1 + 0{zV)) + 2mR , 

Ki^) = ^{W + l)i2W + 1) (1 + Oiz'W')) , 

therefore there exists a solution z^ = iy^ of (f)'_^{z^) = such that 

= -^{W + 1){2W + 1)(1 + 0(i?V(n2iy2))) . 

Under the assumptions of Proposition I3.H all the error terms in (13 ■2p are 
o(l). Also, 

1 + O (exp {-27iNy+}) = 1 + o(l) 



1 — exp {27riA^z} 
uniformly in —1/2 < x < 1/2, since 

These observations and Lemma 1331 justify the saddle point approximation, 
which yields: 



h = {i + o{i)) [— (iy + i)(2iy + i) 
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1/2 f 3K 
exp 



n{W + l){2W + l) 



Now consider the second integral: 

exp q)^\z) p- — — — j -dz 

liy^-l/2 Jiy— 1/2 CXp pTTziVz] - 1 

where = nln/(2;) + 2TTi{R — A^)^. According to Lemma 

•1/2 



I /_ I < exp 0_ ) / exp ^^—cnW"^^ dx 



c 

< ^==exp0_(2?/_ 
VnW 



Choosing ?/_ = — ?/+, we see that 



1 r 3i?2 

|-/_| ^ ^ exp 



v/nVT { n{W + l){2W + l) 



Proof of Proposition \3.S[ According to (13. ip , 

^ ' -N/2<k<N/2 

We have: /(O) = 1, whereas by Lemma [3.51 

I [f{k/N)]'"\ < exp {-cnW^k^N^} . 



Therefore 



^[/(^/Ar)]%xp^ 

|A;|>0 



<^exp<i-^^[^«l 



□ 



□ 



Proof of Proposition \3.3\ (sketch). Consider two cases: ^/nW > N, and y/nW < 
N. In the former case, 

^<j^ E exp{-c*'V/iV^}<C/iV. 

^ ^ -N/2<k<N/2 



In the latter case, we proceed as in the proof of Proposition 13. 1[ making sure 
that one can choose y± so that 

cR^ 



□ 
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4 Embeddings into S 

In this section, we study the number of "topological embeddings" of a graph 
G into S, satisfying certain restrictions. Let us start with the precise defini- 
tions. 

Let G = (V, E) be a connected multigraph, and let it = {ni, ■ ■ ■ , ra^;) be an 
-E-tuple of number^, rig > —1. Construct a new graph G as follows: 

• if ne > 0, replace e = {u, v) with a chain {u, ui, U2,--- , Un^,v); 

• if rig = —1, contract e = {u, v) into a single vertex. 

Finally, consider a system of linear equations C in the variables ni, ■ ■ ■ , n^;: 
(B : ^ Cj(e)ne = n{j) , j = 1, ■ ■ ■ , , 

and denote by Emb(G', (£) the number of subgraphs G" S that are isometric 
to G^ for some it satisfying the system and by Emb+(G, - the number 
of such subgraphs with the additional requirement rie > 0, e & E. 
We impose the following restrictions on (£: 

1. (£ is linearly independent; 

2. Cj{e) >0, I < j < k, e e E; 

3- E'=ic,(e) = 2, eeE. 

These restrictions are satisfied for the systems of equations that we need, 
and slightly simplify the notation. 

Proposition 4.1. In the setting above and under the assumptions 1.-3., if 
< n(l) < n{2) < ■ ■ ■ < n{k) and ^ n{j) = 2n<^W, then 

= (1 + o(l)) N-^-^^ goHl), . . . , nm , 
where go is a continuous homogeneous function of degree E — k, 

go{n{l),---Mk))<\^^. 

The same is true for Emb+(G', €). //X] ^(j) = 2n + 1, 

Emb(G, = Emb+(G', = . 



^to simplify the notation, we write V — #V, E = throughout this section 
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Proposition 4.2. In the setting above and under the assumptions 1.-3., if 
< n{l) < n{2) < ■ ■ ■ < n{k) < min {N^/W\ W), ^ n{j) = In, then 



{2W -\Y 



:i + 0(1)) N (2H^)-^+^-i/ic,(n(l), ■ ■ ■ , n{k)) 



where ho is a homogeneous function of degree M+Y-^JlA^ continuous outside 
the origin, and 

E+V-2k-l 

hG{n{l),--- ,n{k))<^ 



The same is true for Emb+(G', (S). IfYl ^U) = 2n + 1, 

Emb(G, (B) = Emb+(G, (B) = . 
Proposition 4.3. Without any restrictions on n, W , 



i2Wr - iE-k)\^^> +C2Nj^^y^^W 

To apply the results of Section [3], we^ need a simple observation. Let 
u,v e Z/NZ, A,B C Z/NZ. Denote by W„(m, v, A, B) the number of paths 
Pn = UU1U2 ■ ■ ■ Un~iV from u to f in S, such that 

• n 7^ M2, ui 7^ M3, . . . , Un-2 7^ V {''non-backtrackin(f^ ) , 

• Ui ^ A, Un-l ^ B. 

Lemma 4.4. In the notation above, 

f n + 4^A + 4^B\ 

Y^r,{\u-vU)>^n{u,V,A,B)> il-C \/ 1 ^n{\u-vU) • 

Proof of Proposition \4.1\ Choose N'^/W'^ <^ tt-o ^ n{l)/E'^. Let Ag be the 
set of real non- negative solutions of obviously, Ag is an (£"— /c)-dimensional 
polytope with E + k faces. It is not hard to see that 

Vo\E-k-idAe = 0(4^ Vol Ag) . 

nil) 
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Hence, the number n(S) of integer solutions it of ^ with Ue > —1 satisfies: 

n(C) =7(£Vob_,_iA^ (1 + 0(^7^(1))} = Vol A'g{l + 0(EVr2(l))} , 

where Ag = Proj^B-fe-i A^;^ is the projection of Ag onto a set of — /c — 1 
independent coordinates. Moreover, the number of solutions with at least 
one coordinate < uq is at most 

VolAg ■ 0(no^Vn(l)) . (4.1) 

Let it he a. solution with 

ni, - ■ ■ ,nE > no . (4.2) 

It corresponds to 

Ar^-^(2iy-i)-(i + o(l)) 

different embeddings. Indeed, there are iV^ ways to embed the vertices of 
G, and, according to Proposition 13.21 and Lemma [4. 4[ 

(21^-l)"n-^(l + o(l)) 

ways to embed the edges. Thus the total number of embeddings correspond- 
ing to solutions that satisfy fl4.2p is 

Vol A'g N^-^{2W - 1)"(1 + o(l)) . (4.3) 

Applying fl4.ip . Proposition 13.31 and Lemma [4.4[ we see that the number of 
embeddings that violate f l4.2p is negligible with respect to (14. 3p . Therefore 

Emb(G,e) , ,rV-E 



{2w-iy 



(1 + o(l)) VolA'gA^^ 



Finally, observe that Vol Ag is an (ii^ — /c)- homogeneous function of n(l), 
n{k), and that 



VolA'£<Vol{0<Mi,-- - ,UE,Ui + --- + UE<n}: _ 



□ 
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Proof of Proposition \4.2\ First, number the edges 1, ■ ■ ■ ,E so that 1, ■ ■ ■ , V— 
1 is a spanning tree. For e = {ie,te), let Re = i'^ — t'^, where G Ij/NIi 
are the images of ig, te in G' . Note that all the Re are linear combinations of 
■ ■ ■ , Ry-i- To every e G -E we correspond a vector Ve G R^"^ such that 



v-i 

Re = Y.<f)Rf 

/=1 



By Proposition 13.11 and Lemma 14.4 
Emb(G, 



{2W -I] 



,^,,Emb(G,£) 



-Ri,-,-Ri/-i TZ" eeS 



3 



-1/2 



exp 



3 -R^ 



, (4.4) 



where the sum is over rie satisfying Here we have disregarded the contri- 
bution of rie < uq, where 1 -C rio <^ ^(1); this can be easily justified using 
Proposition 13.31 and the asymptotic estimates below. 

Now, applying the Poisson summation formula and discarding the negligible 
terms, we have: 



exp 

iii,-,-Rv-ieZ 



3 Ri 



;i + o(l)) 





^ r 


<^ det 


3 





Ve (» Ve 



-1/2 



According to the Cauchy-Binet formula, 

det ^ CeVe ®Ve= ^ J]^ Ce det^(f e)eer • 



If T is not a spanning tree of G, (fe)eeT is not of full rank and hence has 
zero determinant. If T is a spanning tree, (fe)eeT and its inverse are integer 
matrices, hence the squared determinant is equal to 1, therefore 



det 5^ 



CeVe (» Ve 



T e&T 
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where now the sum is over spanning trees T C E. Going back to fl4.4l) . we 
deduce: 



Emb(G, g) 

(2W^_ l)n 



277 



rin 

It 



-1/2 



rip 



• (4.5) 



Finally, replace the sum ^* with an integral over the polytope 

A'g = ProjRB-fc , 

where Ag is the set of non-negative real solutions of (B, and the projection 
is onto an independent subset of coordinates. This can be done for example 
using the Poisson summation formula. We deduce: 



Emb(G, £) 
{2W-1)" 



2'kW 
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X 



-1/2 



dui---duE^k (4.6) 



T eST 



(where all the other Ue are determined from C) It is easy to see that the 
integral in the right-hand side is ^ -homogeneous in n(l), • • ■ ,n{k). 

Replacing E-p with a single spanning tree Tq, we deduce the upper bound. □ 

Example 4.5. The reader may find the following example illustrative: for 
/c = 1, 

(E: ^ 2rae = 2r 

we have: 



in 



E+V-3 

n 2 



A' 



where the set Ag = n ^Ag = {Y^Xe < IjXg > 0} does not actually depend 
of n. 



We omit the proof of Proposition 14. 3[ which is very similar to the proofs 
of the previous two propositions. 
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5 Proofs of Propositions 12.61 12.7 



According to Lemma [2.51 and the discussion preceding and following it, 
T{n{l), ■■■ , n{k)) = ^ Paths(n(l), ■ ■ ■ , n{ky, ©) , 

where the sum is over connected diagrams ® = (G = {V, E),p), and 

Paths(n(l), ■ ■ ■ ,n{k); d) 

is the number of fc-paths Pn{i),- ,n{k) corresponding to £> and satisfying (a)-(d) 
from Section |2l Let 

, , I times that the ?'-th part I , , 
Cjie) = 4f< r ' >G{0,1,2}, l<3<k,eeE, 

I or p passes through e I 

and consider the system of equations 

: ^ Cj(e)ng = n{j) , l<j<k. 

e 

Then 

Emb+(G', e„) < Paths(n(l),-- - ,n(A;);©) < Emb(G, e^,) ; (5.1) 

here j^V = 2s and j^E = 3s — k. Therefore we can apply the results of 
Section H 

Proof of Proposition \2.(k Choose no and sq so that 

.Tin soAr2 iV2 

If no < n(l) < ■ ■ ■ < n{k) < RN^/^, and £> is a connected diagram with 2s 
vertices and 3s — k edges, 1 < s < So, then Proposition 14.11 yields: 

Paths(n(l), ■ ■ ■ ,n{k); ©) 



{l + o{l))g„{n{l)r-- ,n{k))N-^+' 

(1 + 0(1)) n(l) ■ ■ ■n(A;) fi^,(n(l)/iVV3, . . . ,nik)/N'/') 
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where g^, is homogeneous of degree 3s — k, and 

gt,{xi, ■ ■ ■ ,Xk) = {xi- ■ ■Xfc)"^5fj,(xi, ■■■ ,Xk) . 

If s > So, Proposition 14.31 yields: 

Paths(n(l),---,n(fc);g)) jC'nf^-^'^ (r'N)-^^^ 
(2iy_l)n - {3s-2k)\ ^ ' 

hence the total contribution of all diagrams with s > Sq is at most 

t3 / Ar\s-fc 



s>so 



As the same bound holds for the sum of 

n(l) ■ ■ ■ n{k) {n{l)/N^/\ ■■■ , n{k)/N^'^) 
over diagrams with s > sq, we deduce: 
r(n(l),---,n(A:)) 

where the sum is now over all connected fc-diagrams £> . To prove the bound 
for n{l) < riQ, apply Proposition 14.11 in a similar fashion. □ 

Proof of Proposition \2.'7\ (Sketch). Choose Sq, Uq, Uq so that 

1< < no < W^/^ < < < — , 

and proceed as in the previous proof, using Proposition 14. 21 instead of Propo- 
sition lO □ 
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6 Digression 

Let {fiN}N=i be a sequence of probability measures on R, and let sn ^ +00 
be a sequence of (real) positive numbers. We shall study the scaled measures 

a.,.(A) = 4(l-Ml-2s-A)) 
o-L,^(A) = 4(/^^(-l + 2s^ A)) 

This scaling is meaningful if fiN are close (in some sense) to the Wigner 

measure /xwigner- 

Let 

fT^{n) = / Un{a)dfiN{a) , (6.2) 



where as before Un are the Chebyshev polynomials of the second kind, and 
assume that 



71 

/i7v(ra) = — 



Mn/sN) + i-irMn/sN) +^4^, (6.3) 



'AT 



where 

0i, 0,? e CiO, +00) fl Li (exp { -x-^+' } dx) (6.4) 

(for some 5 > 0), and e]^{n) are "small" (in a sense made precise below), 
for n = 0{sj\f). Then a^^N, o'l,n converge to limits that can be expressed in 
terms of (f)R,(f)L. We state this as a proposition; the main ingredient of the 
proof is a variant of Levitan's uniqueness theorem [17]. The assumptions can 
be definitely relaxed: thus, the integrability condition (16. 4p can be replaced 
with a weaker one using the methods of Levitan and Meiman [18] and Vul 



Proposition 6.1. Let {hn} be a sequence of measures on R, and assume 
that the coefficients fiN (defined in lid. ^) ) satisfy ^6. with (pL, (pR as in 
and 

2. for any i? G IN, eAr(n) = o(l) for n G {1, 2, ■ ■ ■ , Rsn}, with the implicit 
constant depending only on R. 

Then the measures ctr^at and aL^N (defined in lid. converge to limiting 
measures an and (Tl (respectively), where an and ai are uniquely defined 
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by 0/2 and (pi (respectively). The limiting measures a = crR,aL share the 
following properties: 



a(A) - AAf = 0(A) , 



a{X) = O exp<^ -C'|A| — 



l~S/2 



A — i- +00 ; 
A — —00 . 



(6.5) 



Remark 6.2. One can show that 
"+°° sinx^A 



Loo ^^^TT'^'"''^^^ ^ '^''^"'^ = ^^^^^ + 1^^+' 



37r 



3/2 

+ 



(6.6) 



-00 xa/A 

We do not use this in the sequel, and therefore omit the proof. 

To prove Proposition 16. H we shall need the following Erdos-Turan type 
inequality: 

Proposition 6.3 ( [9, Proposition 5]). Let fi be a probability measure on R. 
Then, for any s > 1 and any a G R, 



|/i(a) - yUwigner(a)| < C 

where p{a; s) = max(l — |a|, s^^) 
Proof of Proposition \6.1[ 



+ vp(";s) 2^ 



n=l 



n 



(a) According to assumption 1. and Proposition 16. 3[ 

|rij,7v(A)| <Cmax(A,l) , (6.7) 

where 

2 

tr,n{>^) = o-ij,Af(A) / Vl- a'^da . 

Therefore the sequence {ctr^n} is precompact. The same is of course true for 
{(Tl,n}- If ctr^n, Or, then tr^n^ tr, defined by 

^«(A) = TRiX) + . 
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Therefore we need to prove that a^, are uniquely determined, 
(b) Let 

4n 

(6-8) 

Un+i{a)Ul{a) = c-kUk{a) . 
The exphcit expressions for c^ j^ can be easily derived from the identity 

min(fc,<?) 

Ukia)Ue{a) = ^ f/||._fc|+2m(a) • 

m=0 

We shall only need the following simple properties: 

cl, = 21,^{k) icik/n) + oil))n' , 

4,0 = (co + o(l))n , 

c-fc = 2l2z+i(A:) icik/n) + oil))n^ , 

where c G C[0,4.01] and the o(l) terms are uniform. Thus we have: 
r-4 

xc{x) [(j)R{xn/sN) + 0L(xn/sAr)] dx + CqsI^/tC" + 0(1) , 



n 



f / Un+i{a)Ul{a)d^iN{a) 



4 



xc(x) [0/{(xn/sAr) — (pL^xn/ sjs)] dx + 0(1) . 

Passing to the hmit along a subsequence Nj — )• 00 and applying (16. 7p and 
the dominated convergence theorem, we deduce: 

sin^ x^^J"^ 

- — d(T^(A) = / yc{y)(t>R{xy)dy + cqx''-^ , 

(6.9) 



Ar^ ^^^(^) = / yciy)'pL{xy)dy + cqX 
{xvX) Jo 



-3 
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(c) The relations fl6.9p determine ctr, uniquely. Indeed, according to fl6.7p 
and fl6.4p . the measure must satisfy 



1 + A2 



< +00 









/ exp 






1—00 







daniX) < Ci exp [Cix^-^] ; 



(6.10) 

hence one may apply the argument of Levitan [17]. The latter works as 
follows. Suppose aji is another measure for which fl6.9p also holds. Then aji 
also satisfies fl6.1Up ; = ctr — o^r satisfies 



sm X 



[X 



da{X) = , X e R . 



Set 



...4.y^ 



sm XA 



A2 



dv{X) 



sin^x^A 



A2 



dv{X) 



and apply the Phragmen-Lindelof principle to f{z) in every quadrant. We 
see that / is bounded, hence constant. Therefore Ijr = aR. 

We have proved that the limiting measure is unique; in particular, the 
sequence {aR, N} converges (to this limit.) 

(d) Returning to (lOOD . 

f exp jxv^} daR{X) < Ci exp [Cix^"^] , 

xJ\X\ + Cix^-^^ 



hence 



Taking x — 
follows from ( 16. 7p . 



o-r{X) < Ciexp 



A < ,x > . 



we obtain the second part of (16.51) . The first part 

□ 



Proposition 16.11 can be extended to measures on (for any fixed i.) 
Namely, let {/zat} be a sequence of measures on R^. For simplicity, we assume 
that fiN are symmetric (=invariant under permutation of coordinates). For 

k < i, set 

k 



fj^{n{l), ■ ■ ■ ,n{k)) 



Y\_Un{j)iaj)diJ,N{a) 
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Proposition 6.4. Assume that, for 1 < k < i, 



/i^(n(l), ■ ■ ■ ,n{k)) 



5^ 



IC{l,-,k} 

£jv(n(l), • • ■ ,n{k)) 



+ 



where 



^R,k,^L,k G C((0,+oo)'=) f|i:i(exp [-\\xf-'] dx) 



the coefficients En tend to zero uniformly on it & {1, ■ ■ ■ , Rsn}'', and 

eN{n{l), ■ ■ ■ ,n{k)) 



E 



l<n(l)<---<n(k)<SM 

Then the scaled measures 



n(l) ■ ■ ■ n{k) 



< C . 



(^r,n{\i, ■ ■ ■ ,Xe) = s% [l - fiN{l - 2s^^Ai, ■ ■ ■ , 1 - 2s^^\i)) 
c^L,7v(Ai, ■■■ ,\e) = s% /iAf(-l + 2s^^Ai, ■ ■ ■ , -1 + 2sJ^^\i) 



(6.11) 



converge to limiting measures aR, o"l, which are uniquely determined by {4>k,R}, 
{4>k,L} (respectively). Moreover, a = GL.aR satisfy: 



a 



(Ai,---,A,)-(|,rn,tiA5 



,3/2 



OdlAII^' 



1) 



A,)=0 



exp 



l-i5/2 

-c;iiAii^ 



Aj > , ||A|| —7- oo; 
Aj < , ||A|| — 7- oo 



The proof is similar to the one-dimensional case (Proposition I6.ip : we 
omit it. 



7 Proof of the main results 

Proof of Theorem \l.l[ To show that the random counting measures aR{X), 
(Ti(A) converge in distribution to 2tii(— A), it is sufficient to prove the con- 
vergence of the correlation measures 

peA^i,--- ,A^) = eJJ(Tr(Aj), p^,l(Ai,--- , A^) = E JJo-l(Aj) 
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to 

p,(Ai,-- - ,A,) = Ej]2lii(-A,) 

It will be convenient to denote 

/^iv(") = # jeigenvalues of in (-00, «] \ (7-1) 

(this differs slightly from (11. 5p .) 

Let us first consider l=\. According to Proposition HI 



N 

E^(2ra) = 

2n + l 



0i(2n/Ar^/=^) + (-l)2"0i(2n/Ari/3) 



+ 



N 



0i((2n + l)/iV^/3) _^ (_i)2n+i^^((2n + l)/iVi/3) 



where absorb the difference between the matrices H^^^ / {2Wn — 1)" and 
U2n{H]y/ {2^/2WN), and e^r , 01 are introduced to make the notation compat- 
ible with SectionEl The sequence of measures {/in} satisfies the assumptions 
of Proposition 16.11 with = N'^^^, hence 

Pl,R,Pl,L Pi , 

where pi is a measure determined by 0i (and in particular independent of 
Wn-) For Wn = N/2, pi = pi according to the result of Soshnikov [26] ; 
hence the same is true for any Wn ^ A^^*^^. 
The same argument works for ^ > 1. Indeed, 



n Pen 



where the sum is over all partitions 11 of {1, ■ ■ ■ ,^}. For a subset / C 
{I,-- - ,£}, write n ^ / if 

VP G n P n / G {P, 0} . 



28 



Then 

i=i ^ ^ /c{i,-,^} n-</ 

Now apply Proposition 12.61 and then Proposition 16.41 □ 

Remark 7.1. Another (perhaps, slightly simpler) way to prove the conver- 
gence of the correlation measures is to follow the arguments of /K Section 1. 5], 
and then use the uniqueness theorem for Laplace transform instead of the ar- 
guments of Section 0, as in WU^ . 

Proof of Theorem \1.2[ According to Proposition 12.71 with k = 1, 

■,(1) .(1) 



where again ipi , ejy are introduced to make the notation consistent with 

position with = 



Section [61 Now apply Proposition 16.11 with sn = W^^^ , and deduce that 



ai(A) = ri(A) + —X'f , I ^ dri(A) = ^(x) . 



where 



3vr ' sa/A 

Applying Proposition 12.71 with k = 2 and Proposition 16.41 it is not hard to 
see that 

Eo-L.jv ® aL,N - (Eo-L^jv) ® (EaL,7v) — ^ , 

hence also 

D 

O'R^N, 0'L,N > 0"! • 

□ 

Remark 7.2. Staring at the asymptotics of near zero, it seems natural to 
conjecture that 



^i(A) = f Af + J-|^Af + 0(1) , A->+oo. (7.2) 
OTT V o27r^ 
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We have not been able to prove this as stated. Applying Marchenko 's Taube- 
rian theorem fW^ . one can show that ( [y.jjp holds in a weak sense (say, after 
integrating both sides with a compactly supported twice differentiable kernel.) 

Proof of Theorem M.'d First, for n < Wn, Proposition 14.31 yields: 



EtiH 



(2n) 
N 



N 



E 



Paths(2n; ®) 
{2Wm - 1)" 



5s-4 



Rearranging the sum and using Lemma 12.51 we continue: 



N) 



EtTH 



(2n) 
N 



i2WN - 1)^ 



< 



E 



5s-4 



< Cn < exp 



(3s -2)! 



Ari/2 



N 

H — exp 



(^)! 

^^5/3 



N 



{CsY 



w 



6/5 
N 



w 



2/3 



N 



hence also 

Etrf/2„ 



H 



(2v^ 



< Cn < exp 



+ 



N 



W 



6/5 



exp 



TV 



^^5/3 



2/3 



N 



(perhaps, with a different constant C.) Applying the identities 
have: 



we 



EtrC/. 



H 



2n 



<C <nN + n'^ exp 



H — exp 



6/5 



N 



^^5/3 



2/3 



N 



If A^^/^ < PVat, the right hand side is bounded by Cn^ for n = [N-^^^\ . As 

t/2„,(±(l+£)) 



2n 



> cexp{cn^/e) , 



we deduce that \\Hn/ {2^/2WN)\\ — > 1 (and in fact, 

|(||i/^/(2v/m^) - 1) Ar2/3| 
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AT 



is stochastically bounded.) 

If Wn < N^^\ take n = [W^^^ log^/^ A^J . Then 



and hence again \\HN/{2^y2WN)\\ 1- □ 



8 Random phases 

The steps of the proof for the matrices with entries (11.31) are very similar to 
those for fll.2p . We indicate the necessary modifications. 

• Section [21 Lemma 12.51 remains valid verbatim. Condition (d) in Corol- 
lary 12.21 should be replaced with 

# {{i^j) \ul=U, M^+i = v} = # {(i, j) \ul=V, M^+i = u} . 

In Definition 12. 3[ loops are not allowed, and the third condition should 
be also replaced with 

# {(hj) \ ul = u, ul^^ = v] =^ {{ij)\ul = V, M^^i = u} = 1 . 

Thus, the diagrams in the new sense are a subset of diagrams in the old 
sense. In Lemma 12. 5[ s is now always even, and the estimate is valid 
for even s. In Propositions I2.6|2.7I the sums are now over even s (and 
the functions are different than before.) 

• In Section |5l one should only consider the diagrams that are valid 
according to the new definition, and only even values of s. 

• The argument in Section [7] is still valid. In Remark 17. 2[ we would now 
conjecture that 

a2{X) = + 0(1) , A^+oo. (8.1) 

on 
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9 Concluding remarks 



I. To simplify the exposition, we have only considered the simplest random 
variables (11.21) . (11.31) . Assume that the entries of if at above the diagonal are 
independent, and have symmetric distribution with (uniformly) subgaussian 
tails. Applying the methods of [U Part III], one can assume that 

Ei/i = 1 , HuvER a.s. (0 < |m - v\n < W) 

or 

EHl = , E|ii„,|2 = 1 {0<\u-v\n< W) 

instead of (11.21) or (II. 3p (respectively), and prove analogues of Theorems ll.il 
and 11.21 In particular, this extension covers the frequently considered case 
of matrices with Gaussian elements. 

II. The restriction Wn ^ 1 is also an artefact of the proof. For Wn = 0(1), 
Lemma 14.41 is no longer applicable, hence one should take into account the 
difference between non-backtracking and usual random walk. Thus, we need 
an analogue of Propositions I3.1H3.3I for non-backtracking walks. This can 
probably be proved using either a trace formula for the representation of non- 
backtracking random walk as a Markov chain on the space of directed edges 
(see Smilansky [2^), or the connection to Chebyshev polynomials (see [1].) 

III. It would be interesting to obtain a more detailed description of the 
measure cx/j from Theorem 11.21 In particular, a more precise description of 
the left tail would allow to find the limiting distribution of the maximal 
eigenvalue (cf. the remark after Theorem II. 3p : as to the right tail, it would 
be interesting to justify the asymptotics (17. 2p (and perhaps derive the next 
terms in the asymptotic series.) 

IV. Usual (i.e. non-periodic) random band matrices have non-zero elements 
Huv for 

<\u-v\ < Wn ■ 

We expect that the results of this paper hold, perhaps in modified form, for 
these matrices as well. Following Bogachev, Molchanov, and Pastur [4j, we 
note however that even the limiting spectral measure coincides (II. 6p only if 
1 ^ Wjy <^ N or Wn = (1 — o(l))A^. The limiting spectral measure in the 
complementary regimes has been described by Khorunzhiy, Molchanov, and 
Pastur in [20]. 
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V. It would also be interesting to study the crossover regime 

We refer the reader to the works of Johansson ^3] and Bender [3] for the 
description of the crossover regime at the spectral edge for different kinds of 
random matrices. 

VI. The method of this paper can be used to study the eigenvectors of 
Hjsfl {2yj2WM) that correspond to eigenvalues close to the edge, and, in par- 
ticular, their inverse participation ratio 



If Wn > A^Ve^ the inverse participation ratio of eigenvectors corresponding 
to eigenvalues a = 1 + 0{N^'^^^) is, with high probability, of order A^~^. If 
Wn < N^'\ the inverse participation ratio averaged over eigenvalues in a 
window [1 + a/W'^^ , 1 + h/wH^] is, with high probability, of order W^^^^ . 

VII. We remark that Schenker [22] proved a lower bound on the in- 
verse participation ratio of the eigenvalues in the bulk of the spectrum (for a 
slightly different class of band matrices). In the opposite direction, Erdos and 
Knowles [6l[7] recently proved an upper bound iy^3+°(i) for a wide class of 
band matrices; their argument uses in particular the expansion in Chebyshev 
polynomials developed in the current paper. 

VIII. Finally, there is a natural extension of band matrices to higher-dimen- 
sional lattices: the rows and columns of Hn are indexed by elements of 
{TjjN'LY, and Hn{u,v) = unless < \\u — v\\ < Wn- Similar random 
matrices have been also studied in physical and mathematical literature, cf. 
Silvestrov Disertori, Pinson, and Spencer [5]. We hope to consider the 
spectral edges of such matrices in a forthcoming work. 
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